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We discuss the effects of fluctuations of the local density of charged dopants near a first order phase 
transition in electronic systems, that is driven by change of charge carrier density controlled by 
doping level. Using a generalization of the Imry-Ma argument, we find that the first order transition 
is rounded by disorder at or below the lower critical dimension dc = 3, when at least one of the two 
phases has no screening ability. The increase of dc from 2 (as in the random field Ising model) to 3 
is due to the long-range nature of the Coulomb interaction. This result suggests that large clusters 
of both phases will appear near such transitions due to disorder, in both two and three dimensions. 
Possible implications of our results on manganites and underdoped cuprates will be discussed. 
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Study of effects of quenched disorder on first order 
phase transitions has a long history. The prototype 
model for such studies is the random field Ising model 
(RFIM) [1]. At low-temperature the pure Ising model 
spontaneously develops magnetization and there are two 
different phases corresponding to the two possible direc- 
tions of the magnetization. A uniform external field h 
drives the systems through a first order phase transition 
at he = 0, where the magnetization changes discontin- 
uously. If a random field with zero mean is applied on 
the other hand, it would like the local order parameter 
to point in the direction of the local field, thus clusters or 
domains of both phases would form, the global magneti- 
zation would be destroyed, and the first order transition 
associated with the jump of global magnetization would 
be rounded. In order for this to happen however, the 
energy gain due to forming clusters with magnetization 
following the local random field direction. Eh, must be 
large enough to overcome the energy penalty associated 
with the domain walls separating the clusters with oppo- 
site magnetization, E^. Imry and Ma [1] analyzed the 
competition between these two effects as a function of the 
linear size of the cluster L, and found that for systems 
with dimensionality d less than the lower critical dimen- 
sion dc = 2, Eh ~ L'*/^ overwhelms E^ ^ L^^^ for large 
enough L; in this case clusters will form at sufficiently 
large length scales no matter how weak the random field 
is, and the first order transition is rounded. For d > dc 
on the other hand the global magnetization and first or- 
der transition remain intact for weak random field. For 
d = dc = 2, Eh and i?i„ scale with L the same way: 
Eh ^ i?i„ ^ L. More detailed studies [2,3] indicated that 
at d — dc ~ 2, Eh eventually dominates i?i„ at large L 
and destroys the ordered phase and the associated first 
order transition. These considerations and conclusions 
have been generalized to generic first order phase transi- 
tions [4]. 

Recently there has been considerable interest in the 
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FIG. 1. Schematic phase diagram of an electronic system. 
As the doping level x is varied, there are two phases separated 
by a first order phase transition ai x = Xc- 

rounding of first order phase transitions in electronic sys- 
tems, and the associated cluster formation. In these sys- 
tems the first order transition in the absence of disorder 
is driven by the change of charge carrier density, which 
in turn is controlled by the density of dopants introduced 
into the system; see Fig. 1. Here the dopants play a dual 
role. First of all they provide the charge carriers (con- 
duction electrons or holes), and their density determines 
the location of the system in the phase diagram. On the 
other hand they are also the source of disorder, as they 
distort the local lattice structure and more importantly, 
there is local fluctuation in their density. As we argue 
below, due to the long-range nature of the Coulomb in- 
teraction between the dopants and the charge carriers, 
the effect of dopant density fluctuation is much stronger 
than the random field in the random field Ising model; 
the latter only couples to the local order parameter and 
thus represents short-range interaction. A consequence of 
that is it increases the lower critical dimension dc from 2 
to 3; we thus expect in both 2- and 3-dimcnsional systems 
the otherwise first-order phase transitions are rounded by 
dopant density fiuctuations, and clusters of both phases 
will always form near the phase boundary. In the fol- 
lowing we illustrate these points by performing an ex- 
tension of the original Imry-Ma analysis. We will con- 
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sidcr two different generalizations of the Coulomb inter- 
actions in 3D to arbitrary dimensions. The first is to keep 
the momentum dependence: V{q) ~ l/q^, which yields 



in real space (in 2D, V^^''{r) ^ logr): the 



other is simply using V^^\r) ~ 1/r in all dimensions. We 
will discuss both cases. 

Now assume that the average dopant density of the 
system is n = ric, i.e. at the transition point (or doping 
level X = Xc 'vo- Fig. 1). In the absence of any disorder, the 
system will be in one of the two phases depicted in Fig. 
1 and there are no clusters or domains, since generically 
there is an energy penalty associated with the domain 
walls that is proportional to the area of the walls. We 
now discuss what happens when there is local fluctua- 
tion of the dopant density. Consider a domain of linear 
size L. The average number of dopants in this domain 
is N{L) ~ L'^, while the typical fluctuation of N{L) is 
5N{L) ~ y/WV) ~ Obviously, the local fluctu- 

ation of dopant density tend to force domains of both 
phases to form, according to the local density of dopants; 
the energy penalty of not forming domains is dominated 
by the electrostatic energy due to the imbalance between 
the dopant charge and carrier charge in the regions that 
are in the "wrong" phase where the carrier density can- 
not reach that of the local dopant density, which is 



Ec{L) ~ [5N{L)fV{L) ~ L'^V{L) 



(1) 



per domain. For the two types of generalization of 
Coulomb interaction, we find 



and 



(2) 
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respectively. In order for the domain formation to oc- 
cur, the gain of Coulomb energy Ec must overcome the 
domain wall energy penalty 



rd-l 
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Comparing the expression of Eyj with those of Ec^^' and 
, we find for the first generalization of Coulomb in- 
teraction, Ec^^ dominates E^ at sufficiently large L for 
d < dc = 3, i.e., the lower critical dimension above which 
the first order transition remains intact is rfc = 3; at 3D 
the effect of the dopant density fluctuation is marginal. 
For the second generalization of Coulomb interaction, 
Ec scales with L exactly the same way as E^, regard- 
less of the dimensionality d. Thus the effects of dopant 
density fluctuation is always marginal. 

These results are very different from that of the ran- 
dom field Ising model (RFIM), in which case dc = 2; 
the origin of the fundamental difference lies in the long- 
range nature of the Coulomb interaction, which yields 



very different dependence of Ec on the domain size L, as 
compared with Eh ~ L'^^^ of the RFIM. 

The discussions presented above are quite general. In 
the following wc propose a simple lattice model that con- 
tains the basic ingredients of the physics being considered 
here, which can be studied numerically to check the re- 
sults obtained above [5]: 



H 



(pi - ni){pj - Uj) J 



^ [1 - sgn{pi ■ pj] 



<ij> 



(5) 



Here i and j are indices of lattice sites of a d-dimensional 

cubic lattice, Ui is the quenched background (or dopant) 
charge which follows some specific distribution like Gaus- 
sian or box distributions, and pi is the electron charge 
which is the dynamical variable of the model, and it is 
varied to minimize H. The first term represents the 1/r 
Coulomb interaction among the charge carriers, and the 
summation is over all pairs of sites in the system; it can 
also be replaced by other forms of long-range interaction 
that is appropriate. The second term represents domain 
wall energy cost and the sum is over nearest neighbors; 
here without losing generality we have assumed the phase 
boundary is located at pc = 0, thus there is an energy 
cost of J for each nearest neighbor bond that connects 
sites with opposite signs of p, or sites occupied by dif- 
ferent phases (the function sgn returns the sign of its 
argument). It is clear that in the absence of fluctuation 
in the background (or dopant) charge density, i.e., Ui = n 
independent of i, we will simply have Pi = n everywhere, 
there is no competition between the two terms in H, and 
there is a first order transition at n = from the phase 
with /> < to the phase with p > 0. On the other hand 
when Hi follows certain distribution with a finite width, 
and its average n is close to 0, these two terms compete 
with each other. The first term prefers to have pi match 
Ui everywhere, thus pi with both signs (or both phases) 
would appear in the system; the fluctuation of Ui plays 
a role similar to the random field in the RFIM. On the 
other hand the second term penalizes near neighbor pairs 
with opposite signs of p; J plays a role similar to the Ising 
coupling in the RFIM. Numerical simulation of this and 
other closely related models are currently underway [6]. 
A few comments on our results arc now in order. 

(i) It is very important in our discussion that the first 
order transition is driven by the change of doping level. 
There are, of course, other first order transitions in elec- 
tronic systems that are driven by other mechanisms (like 
change of pressure) without changing the doping level; 
our results do not apply to those cases, as long as the 
position of the first order transition does not involve the 
doping level. 

(ii) Within the RFIM, it is known that if there ex- 
ists power-law correlation in the random field, the lower 
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critical dimension dc may be increased [7]. The increase 
of dc due to long-range Coulomb interaction discussed 
here has some resemblance to this effect, but the precise 
mechanism of increase of dc is different. 

(iii) In the case of the RFIM, it has been proven 
rigorously [3] that right at the lower critical dimension 
d = dc = 2, the random field destroys the ordered phase, 
thus the first order transition is rounded by the cluster 
or domain formation [8]. We are unable to provide such 
a proof in the case we study here. However heuristically 
it is easy to see why a,t d= dc disorder effects eventually 
dominate. This is because for a given domain size scale 
L, there is always some probability that the fluctuation of 
the number of dopants within a domain is strong enough 
to overcome the energy penalty of the corresponding do- 
main wall; such probability does not decrease with in- 
creasing L sX d = dc- Therefore domains with arbitrary 
size will appear at d = dc- We thus expect that for the 
1/r Coulomb interaction, disorder effects dominate and 
the first order transitions are rounded by cluster forma- 
tion in both d = 2 and d = 3. 

(iv) Thus far we have not explicitly discussed the na- 
ture of the phases that are involved, and how they may 
affect the consideration above. If both of the compet- 
ing phases are insulating and have no screening ability, 
then all of our discussion will just come through with no 
modification [9]. The situation is more complicated and 
interesting when metallic phases are involved. Metallic 
phases have screening abilities; they can maintain charge 
neutrality of the system at large scales by creating lo- 
cal carrier density deviation near a charged impurity to 
compensate for the its charge; locally^ the carrier den- 
sity may deviate from its average value significantly and 
may even move across the phase boundary. Thus when 
both phases are metallic, the charge imbalance discussed 
above may be compensated for by screening processes, 
and our conclusions do not apply beyond the screening 
length. This is not the case however when one of the 
two competing phases is insulating; this is because the 
insulating phase does not have the screening ability, thus 
cluster formation cannot be avoided this way when the 
phase boundary is approached from the insulating side, 
and the transition from an insulator to a metal is a perco- 
lation transition. Even when both phases have screening 
abilities, our results are still relevant when the screening 
length ^ is long compared to microscopic length scales 
(typical for " poor" metals). In this case clusters will 
form at d = 3 up to size i ~ ^ no matter how weak the 
dopant fluctuation is. 

We emphasize that the considerations presented here 
arc rather generic, and the results should be relevant to 
any first order phase transitions driven by change of dop- 
ing level in electronic systems. As possible examples of 
the situation discussed here, we now turn our discussion 
to two electronic systems that are of tremendous current 
interest, namely the manganites with colossal magneto 



resistance and high Tc cupratcs. Both compounds have 
very rich phase diagrams, and in particular, there is a 
transition from an insulating phase to a metallic phase 
as dopant concentration is varied (in the latter case the 
metallic phase is superconducting at low temperature). 
Recently there has been considerable experimental ev- 
idence suggesting that clusters of both phases co-exist 
and the system becomes inhomogeneous, when the tran- 
sition is approached [10-17]. It has been suggested that 
such inhomogeneity is induced by disorder near an oth- 
erwise first order transition [18-22]. In the case of the 
manganites, convincing arguments have been advanced 
that the cluster formation and the percolative nature of 
the metal- insulator transition may be the origin of the 
colossal magnetoresistance [18-20,22], and prediction of 
possible similar behavior in the cupratcs has been made 
based on the similarity between these two systems [22]. 
We note that the manganites are 3D systems, which is 
above the lower critical dimension of the RFIM, dc = 2. 
The cuprates are quasi-2D systems; however its large dis- 
tance behavior must be determined by the 3D nature of 
the system. Thus without long-range interactions, the 
formation of clusters and the resultant inhomogeneity 
due to disorder near a first order transition is possible 
but would not be guaranteed in either case [23]. The 
importance of the long-range nature of the Coulomb in- 
teraction, as we discussed here, is to guarantee cluster 
formation and inhomogeneity by pushing the lower crit- 
ical dimension up to dc = 3. Our results are thus consis- 
tent with the phenomenology of both the manganites and 
cuprates, and point to the generality of cluster formation 
near first order phase transitions in electronic systems, 
which is of central importance in the theory of Dagotto 
and co-workers on the colossal magnetoresistance in the 
manganites, based on phase separation and cluster for- 
mation. 
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